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The complexity of a signal can be measured by the Recuri'ence period density entropy (RPDE) from the reconstructed phase space. 
We have chosen a window based RPDE method for the classification of signals, as RPDE is an average entropic measure of the 
whole phase space. We have observed the changes in the complexity in cardiac signals of normal healthy person (NHP) and 
congestive heart failure patients (CHEP). The results show that the cardiac dynamics of a healthy subject is more complex and 
random compare to the same for a heart failure patient, whose dynamics is more deterministic. We have constructed a general 
threshold to distinguish the border line between a healthy and a congestive heart failure dynamics. The results may be useful for 
wide range for physiological and biomedical analysis. 
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1. Introduction 

Human heart reveals electrical discharges from specific lo¬ 
calized nodes within the myocardium. These discharges propa¬ 
gate through the cardiac muscles and stimulate contractions in 
a coordinated manner in order to pump deoxygenated blood via 
the lungs (for oxygenation) and back into the vascular system. 
The physical action of human heart is therefore induced by a lo¬ 
cal periodic electrical stimulation. As a result of that, a change 
in potential can be measured during the cardiac cycle by elec¬ 
trodes which are attached to the upper torso of heart (usually 
both side of the heart). This recorded signal is known as elec¬ 
trocardiogram (ECG) [1]. A typical ECG waveform, consists 
of P,Q,R,S,T,U & V major components, reveals definite pattern 
in the oscillation. However, different pattern can be observed 
due to changes of heart condition. The increasing and decreas¬ 
ing mechanical activities are the major causes for this complex 
phenomena in heart dynamics [2]. This paper studied the com¬ 
plexity of ECG signal of NHP and CHEP using the method of 
phase space analysis. 

Phase space analysis is one of the most useful methods for 
explanation of long term dynamics. It is an abstract Euclidean 
space that reflects asymptotic nature of the interconnected vari¬ 
ables which are responsible for the original dynamics [3-5]. 
The number of such variables is known as Embedding dimen¬ 
sion [4] in which trajectory of the phase space can be flour¬ 
ished properly. Eor a continuous signal, reconstruction of phase 
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space can be done by suitable time-delay and proper embed¬ 
ding dimension [6]. Suitable time-delay is generally obtained 
by the method of Average Mutual Information [7,8] and proper 
embedding dimension is obtained by method of Ealse near¬ 
est neighbor [9-11]. However, different types of trajectory’s 
movements have been observed in the phase space, viz; peri¬ 
odic, quasi-periodic, chaotic,etc, which can be described from 
Recurrence plot (RP) [12-16]. 

Recurrence plot is a diagrammatic representation of a 2D ma¬ 
trix, whose elements are considered 1 if it is recurrent; oth¬ 
erwise it is considered 0. 1 and 0 represented by black and 
white dots respectively in RP. RP quantihes the structure of the 
phase space; it does not have anything to do directly with the 
signals. The diagonal lines in RP are used to measure the com¬ 
plexity [13]. But those lines reflects the parallel movements of 
the trajectories, which explain the deterministic nature of the 
system [13]. So measuring complexity using diagonal lines in 
RP characterizes the degree of chaos present in the phase space 
dynamics. Recently, another measure of complexity with a dif¬ 
ferent form of recurrence have been proposed in [17], which 
quantifies the presence of non-deterministic term in stochastic 
dynamics. So far, there is no particularly effective tool to quan¬ 
tify the complexity of a stochastic signal. Therefore, to investi¬ 
gate this type of complexity, we have implemented Recurrence 
period density (RPD) method of reference [18]. 

RPD is a probability measure which calculates density of the 
recurrent times. In RP, recurrent times are considered as length 
of white horizontal/vertical lines between each pair of recur¬ 
rent points. It can be observed that, a system possesses more 
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complex dynamics if length of the white lines varies frequently 
in RP. Consequently, variation of recurrent times becomes very 
high. Thus, probability of recurrent times reveals a proper jus¬ 
tification to identify the complex dynamics. The quantification 
parameters Recurrence period density entropy (RPDE) mea¬ 
sures complexity of such dynamics. In this concern, Shanon 
entropy is utilized by probability density of recurrent times. So 
far, the analysis of complexity was done with respect to the cor¬ 
responding time series only [19-21], which does not predict 
the long-term behavior of the dynamics. So, to investigate the 
long-term behavior, it is important to study complexity of the 
phase space for more accurate prediction. 

It is well known that a deterministic system is a model, which 
always produces the same output for a particular initial state. 
Any ODE system (autonomous/non-autonomous) ^ = FiX, t) 
is deterministic by nature and can produce chaotic scenario with 
proper conditions. The most important point is, there are non- 
deterministic (random/stochastic) systems, generated by dif¬ 
ferential equations, produce chaos. Eor example, the additive 
noise term (^) with a deterministic model ^ = F(X, t) + ^(t) 
can make the whole system non-deterministic [22], as the same 
initial condition does not give same output. In fact, this is 
an example of a non-deterministic system produces stochas¬ 
tic chaos. This kind of chaos can be observed in Ecological 
models. Lasers [23-25] and other semiconductor devices. In 
many cases noise can increase the nature of the complexity of 
the system, sometimes it may be useful to decrease chaos [26] 
to revert back the system to a regular state. Although, a noisy 
signal or a purely stochastic signal does not have any determin¬ 
ism and can be investigated properly using the corresponding 
time series, sometimes the non-deterministic outputs can also 
be generated by the governing equations. Noise induced (also 
added) chaos are very interesting to investigate, due to its rich 
complexity and non-deterministic nature. Recently there are 
research on noise induced synchronization and also chaos syn¬ 
chronization between stochastic models. Even the Chua system 
can be stochastic [27] with induced noise, studied both theoret¬ 
ically and experimental observations. In this article, we have 
proposed the dynamical complexity of some known nonlinear 
deterministic and stochastic systems as well as power noise by 
RPDE analysis. Eurther, by defining window Normalized Re¬ 
currence period density entropy (NRPDE), we have categorized 
ECG of NHP and CHEP into two classes. In fact, mean window 
NRPDE successfully distinguishes both types of ECG in term 
of complexity. Eurther, a proper threshold can be found in the 
mean window NRPDE by which we can conclude that whether 
a NHP becomes CHEP or not. 


2. Recurrence plot and Recurrence period density 

Eor a n-dimensional phase space X = [(x,) : x, e 51”,/ = 
1,2,..., N], recurrence means closeness of any two points. Two 
points are considered close if their state vectors certainly lie in a 
e-neighborhood. Eormally, two points x ,,xj e X,i - 1,2 ,...,N 
are recurrent if ||x, - Xj\\ < e. The recurrent matrix is thus 


defined as 


Ry = 0(e-||xi - Xjll),/= 1,2,...,A, (1) 

where © is the Heaviside function, ||.|| is Euclidean norm of 
the reconstructed phase space, and e is radius of the neighbor¬ 
hood. RP corresponds recurrent and non-recurrent point by ‘1’ 
(black dots) and ‘0’ (white dots) respectively. Phase space is 
reconstructed by suitable time-delay (t) and proper embedding 
dimension (m). An m-dimensional reconstructed is given by 


((« 


( 2 ) 


Appropriate choice of embedding delay t plays an impor¬ 
tant role in this context. Selecting t too small means that any 
trajectory lies close to the diagonal of RP and hence spurious 
recurrences appear. Eor too large t, state vectors in the em¬ 
bedded space fill a large cloud. In that case, recurrences will 
be difficult to find out. The optimum time-delay is chosen by 
AMI method [7,8]. AMI method measures average information 
shared by signal itself. Eor a signal x(/), /, 1,2,.., N and time- 
delay/lag (t), ami is defined by 


4 \ - V x{i + t)) 

^ Pr{x{i))Pr{x{i + t)) ’ 


( 3 ) 


where Pr{.) denotes probability. In order to select suitable r, 
we used the first local minimum principle which is proposed by 
Eraser and Swinney in [7]. 

Embedding dimension reveals number of necessary indepen¬ 
dent coordinates for phase space reconstruction. In order to 
estimate embedding dimension, we used the method of false 
nearest neighbours (ENN). This method is based on the idea 
of Ealse neighbour state. Ealse neighbour state occurs for self 
crossing trajectory in a phase space. In fact, when a trajec¬ 
tory is projected to a space with too small dimension, trajec¬ 
tory crosses itself and the so called false neighbour states oc¬ 
cur. As dimension of phase space increases, number of trajec¬ 
tory self-crossings and false neighbours decreases. It has been 
observed that, both disappear completely for too large dimen¬ 
sion of the phase space. To determine ENN for a time series 
x(/), i - 1,2,.., N, we used the following criterion; 

\x(i + rriT) - x^^{i + mr)! 

1-2 - 11 - 2 > A, ( 4 ) 

R 

where x^^ represents points of nearest neighbourhood (NN). 
R is given by R^ - < x >]^, where < x >= 

jj x(i)- As per [9], A is so chosen that it stays around 
A = 2. 

Now RP reveals various structure that provides information 
about the nature of phase space [13]. Parallel movements 
and trapping situation/laminar states are described by diagonal 
lines and vertical/horizontal lines respectively. Periodic/quasi- 
periodic phase space nature can be understood by the presence 
of only diagonal lines with equal/unequal time span. Rectan¬ 
gular like structure consists of diagonal line with some isolated 
points and vertical lines represents chaotic regime. However, 
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more complex structure can be observed in RP for stochastic 
chaotic or purely stochastic system. Uniformly distributed iso¬ 
lated points corresponds uniformly distributed white noise. Sin¬ 
gle isolated points imply strong fluctuation exists in the process, 
i.e; the process is either uncorrelated or anti-correlated. Previ¬ 
ous characteristics can also be verified by Recurrence period 
density. 

The measure of RPDE is developed from the idea of recurrent 
time between the recurrent points [18]. The Recurrent times 
are calculated by counting non-recurrent points or white lines 
between two recurrent points Xi,Xj in the RP Rij. Thus, for any 
two recurrent points x,, xj e Rij, the recurrent time denoted by 
Tt, is defined as = (/ - j). In fact, Ty corresponds to the 
least recurrent time,72 corresponds to the next and so on. Now, 
for all points in Rij a series of recurrent time interval R{Tk) 
is obtained as the number of occurrence of Tt- RPD P{Tk) is 
finally defined as the probability of R{Tk) among the sample 
space {R{Ti^}. This is given by (5). 


PiTk) = 


R(Tk) 


( 5 ) 


where T^ax - max[Tk}. 

The changes of RPD as observed from RP are numerical il¬ 
lustrated for some commendable complex systems- 
Lorenz system. Logistic map and Stochastic system- first or¬ 
der ordinary differential equation with additive noise. We com¬ 
puted PiTk) for each system by (3) and observed their individ¬ 
ual characteristics. This analysis is also done for one dimen¬ 
sional colored noise which is generated by ^-law. The whole 
computational result reveals a clear conception about the pat¬ 
tern of probability density which changes due to the increase of 
complexity of the system. 

Lor the Lorenz system- 


dx 

— =a{y- x), 
at 

(6a) 

1 

II 

^1^ 

(6b) 

dz 


(6c) 


(a - 10, r = 28,x(0) = 8,y(0) = 9,z(0) = 25) probability den¬ 
sity is calculated for two different conditions: double periodic 
(b - and chaotic (b = |). Since most of the state vec¬ 
tors lie in parallel trajectories for b - variation of recurrent 
points are very few in number (Lig.la). Consequently, proba¬ 
bility density appears for few recurrent times (Lig.lb). Lorenz 
system ensembles chaotic regime with heteroclinic transition 
for b - |. As a result, lesser state vectors can be found in 
same parallel trajectory and different types of recurrent time are 
found in the corresponding RP (Lig.lc). Corresponding proba¬ 
bility density is shown in Pig.Id. 

In logistic system: 


x„+i = ax„(l-x„), (7) 

we consider a = 3,4. The concept of parallel movement of 
trajectory is completely meaningless here. In RP, black dots 


represents only recurrence in e-neighborhood. Lor a - 3 and 
4, corresponding RP and probability distribution are given in 
Pig.le, Ig, If and Ih. It is seen that, the probability of recurrent 
time occurs in a very short range of recurrent times (Pig. If) for 
a - 3. On the other hand, various RPD appears for large range 
of recurrent time for a - A (Pig.lh). Therefore, the system (7) 
is less deterministic when a - A. 

Next consider a Stochastic differential equation (SDE) with 
exponential decay: 


^ = -Zr H- cr(z,)^(f), (8) 

at 

where ^(f) is additive random noise and crizt) is strength of the 
noise. 

Lor crizt) - 0.5 and 1.5, (T,m) are found to be (17,5) and 
(15,4) respectively. The time-delay r and embedding dimen¬ 
sion m are calculated using (3) and (4) respectively. Lor the 
purpose of unavoidable restriction the corresponding RP are not 
given in this manuscript. 

We have also investigated the natures of RP and RPD from 
^-noise via phase space reconstruction (2). Random or dis¬ 
ordered structure is observed in RP of the noisy signal for 
y6= 0(Pig.2a). In this case, the length of variation of recur¬ 
rent times is higher than that of the same obtained for other 
values of /?. On the other hand, an well patterned RP is ob¬ 
served in Pig.2d for fi- 2, which implies low recurrent time 
variation. The same investigation has also been done for the 
intermediate cases- /?= 0 (Pig.2b),y6= 1.5 (Pig.2c), which indi¬ 
cates the gradually well structured RP. Thus, study of RP on 
^-noise suggested that the variation of recurrent times are in¬ 
versely proportional to the power of the noise fi. The nature 
of probability density can be observed from each RPD for ^- 
noise (J3 - 0, 0.5,1.5,2) with Tk,k — 1,2,.., 1000 (see Pig.2e). 
It is seen that, values of PiTk) follows decreasing trend with 
Tk,k e [1,1000] for ^-noise ij3 - 0,0.5,1.5,2). In fact, the 
tendency of getting lower values for PiTk) increases as fi in¬ 
creases (e.g; see the values of PiTk) for green and yellow). It 
implies that the variation of recurrent times, i.e; Tk decreases as 
P increases. The same behaviour is observed for intermediate 
cases-yS = 0.25,0.75,1,1.25,1.75. This indicates that complex¬ 
ity is inversely proportional to the power noise. 

RPD quantifies the complexity of the phase space with re¬ 
spect to the variation of parameters. But RPD can not mea¬ 
sure the order of complexity. The complexity of a system 
can be quantified by Normalized Recurrence period density en¬ 
tropy (NRPDE). We have extended the idea of NRPDE to win¬ 
dow NRPDE by taking some partitions of the solution com¬ 
ponents/signals to observe the variation of complexity within 
different time frames. 

2.7. Window Recurrence period density entropy 

Recurrence periodic entropy (RPDE) is based on the concep¬ 
tion of Shanon entropy, which is defined as 

N 

^7 = - ^ Pixdlogipixi)), (9) 

/=1 


3 




(a) 



(b) 
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Figure 1; (Color online) (a) represents the RP of the system (6) for chaotic phase space (with e = JqD, D is the diameter of the 
phase space). The RPD for both = | are shown in (b) (with 1000 recurrent times) respectively, (c) represents RP of the system 
(7) for a- 3 (e = 0.1 cr, cr is the standard deviation of the solution). The corresponding RPD are given in (d) (for 3). 
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Figure 2: (Color online) (a)-(d) represents RP of ^-noise, for fi- 0,0.5,1.5,2 respectively.(e) represents log - log plot of RPD of 
^-noise, fory8= 0 (green), 0.25, 0.50 (red), 1.5 (black), 2 (yellow). Using (2) and (3), reconstruction parameters (t, m) are found 
to be (3,4) (forjS = 0), (5,4) (for^S = 0.5), (14,5) (for^S = 1.5) and (31,5) (for/I = 2). In each case, thresholds e are chosen as 
e - O.lcr, where cr is the standard deviation of the power noise. 
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where p(xi) is the probability of the event x, (by convention 
QlogO is taken as 0.). Hence, RPDE is given by 

T„ax 

// = -£p(r,)iogP(r,). (10) 

k=\ 

Since T^ax changes with sampling time, so a normalization is 
necessary for RPDE. Thus, normalized RPDE is defined as 

T,r,ax 

finorm ~ -ilogT^axT' 2 P(Tk) logP(r,), (11) 

k=l 

where log P(i) ^og P(i). 

H„orm actually RPDE of the reconstructed phase space, where 
the points are independently identically distributed. 

Eor Window NRPDE, consider a time series {x(/))j^j and the 
intervals //, = [x(l + (A^ - \)L), x(l + NL)], 7 e Z. Eor each //,, 
reconstruction of attractors are done by (1). Let ti and nii are 
the embedding parameters for (1). Then, using (9), we define 

tLx 

HLm = -ilogTk,aX^ £ P(T^)logP(T^), (12) 

7’i=l 

where and denotes the recurrent times and maximum 
of them over the interval Ii. describes the normalized en¬ 
tropy of the phase spaces which are reconstructed over window 
Ii. Window NRPDE measures variation of complexity 

of the dynamics over different time frame. 

NRPDE is calculated for different values of the parameter 
for the Lorenz system, logistic map and first order SDE. The 
changes of the NRPD values is given in Eig.3. It has been ob¬ 
served that, with the increasing values of the parameter, the cor¬ 
responding NRPDE also increases (Eig.3a). Similar nature can 
be observed for logistic map (Eig.3b). Eor a€ [3,3.5], the val¬ 
ues of NRPDE is either zero or tends to zero as evident from 
Eig3b . In this range, the logistic map executes periodic dy¬ 
namics. Also values of NRPDE increases as a increases. So, 
changes of the nature of complexity of the phase space for 
Lorenz system and Logistic map can be measured by NRPD 
with parameters. In case of stochastic differential equations, 
decreasing patterns of NRPD can be observed with increment 
of the strength of the random noise(Pig.3c). It signifies that 
whenever strength of the noise increases, dynamics become less 
complex. Thus, NRPDE is an effective tool to describe the na¬ 
ture of complexity of the dynamics. 

However, changes of NRPDE with parameters can not be cal¬ 
culated for real data. To observe the changes of dynamics com¬ 
plexity of a given signal, it is effective to investigate the NRPD 
in different time frame. 

3. Recurrence plot and Recurrence period density 

5.7. Application of window NRPDE on ECG signals 
3.1.1. Collection of Cardiac signals 

Two types of ECG signals-(fl) Normal healthy persons 
(NHP), (b) Congestive Heart Eailure Patients (CHEP) are con¬ 
sidered as experimental subjects. All ECG signals are down¬ 
loaded from Combined measurement of ECG, breathing and 


seismocardiogram (CEBS) database and BIDMC Congestive 
Heart Eailure Database available in Physionet [28]. 

The CEBS database contains 20 presumed healthy volunteers 
(12 men, aged 19 to 30, and 8 women, aged 22 to 28) who were 
asked to be very still in supine position on a comfortable sin¬ 
gle bed and awake. This subjects are having healthy (8 men, 4 
women) and sedentary (4 men, 4 women) life style. It is also 
observed that, all volunteers are non-smoker. Signals are ac¬ 
quired by Biopac MP36 data acquisition system (Santa Barbara, 
CA, USA) with a bandwidth [0.0577z, 150/7z]. Each channel 
was sampled at 5 kHz. 

The BIDMC Congestive Heart Eailure Database contains 15 
long-term ECG signals (11 men, aged 22 to 71, and 4 women, 
aged 54 to 63) with severe congestive heart failure (NYHA class 
3-4). This group of subjects was part of a larger study group re¬ 
ceiving conventional medical therapy prior to receiving the oral 
inotropic agent, known as milrinone. The individual record¬ 
ings are each about 20 hours in duration, and contain two ECG 
signals each sampled at 250 samples per second with 12-bit 
resolution over a range of 10 millivolts. The original analog 
recordings were made at Boston’s Beth Israel Hospital (now 
the Beth Israel Deaconess Medical Center) using ambulatory 
ECG recorders with a typical recording bandwidth of approxi¬ 
mately 0.1 Hz to 40 Hz. Annotation files (with the suffix .ecg) 
were prepared using an automated detector and have not been 
corrected manually. 

3.1.2. Complexity analysis 

Complexity of NHP and CHEP can be analyzed from its cor¬ 
responding reconstructed phase spaces. As a primary task, RPD 
is calculated on whole signals for testing complexity of the dy¬ 
namics. During the reconstruction (2) of ECG signals, it is 
observed that most of the time-delay t e [20,30] (for NHP), 
T e [12,19] (for CHEP) and embedding dimension m e [3,6] 
(NHP), m e [5,9] (for CHEP) respectively. As a sample illus¬ 
tration, RP of the ECG signal of one NHP and one CHEP and 
RPD of all NHP’s and CHEP’s are given in Pig.4. RP of NHP 
(Pig.4a) shows more bowed lines than diagonal lines. Long 
bowed line structure indicates that the evolution of states is sim¬ 
ilar at different epochs but with different velocity, which results 
in rapid change of the dynamics [13]. However, these bowed 
lines are too short and also very few recurrent points appear in 
the RP of NHP. This implies that the cardiac dynamics is less 
deterministic and hence more complex for NHP. On the other 
hand, RP of the CHEP (Pig.4b) shows many diagonal lines, also 
many recurrent points, which reveals that the cardiac dynamics 
is more deterministic and hence less complex. These can be 
best judged from RPD analysis. It can be observed from the 
RPD analysis of NHP that probability density P{T) appears for 
most of the time recurrences T (Pig.4c). Thus, variable recur¬ 
rent times can be observed in RP. On the other hand, PiT) exists 
for less recurrent times (Pig.4d) for CHEP. So, recurrent times 
is low in the phase space and hence in the corresponding RP. 
Therefore, movements of the trajectories in the phase space of 
NHP are more complex than that of CHEP. The corresponding 
window NRPDE is given in Eig.5. In Eig.5a, window NRPDE 
for NHP shows always higher values (see color bar). The higher 
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Figure 3: (Color online) NRPD of (a) system (6) with r e [11,70], (b) system (7) with ae [3.001,4], (c) system (8) with cr(z,)e 
[0.1,1.5], 



Figure 4: (Color online) (a) and (b) represents RP of one NHP (with (r, m) - (20,4)) and one CFIFP signals (with (r, m) - (12,6)) 
respectively. In each case, thresholds e are chosen as e = O.lcr, where cr is the corresponding standard deviation of the ECG signals. 
RPD of (a) NHP with 1000 recurrent time, (b) CHFP with 1000 recurrent times. For NHP embedding dimension and time-delay 
are calculated by FNN-method and AMI-method respectively. Different colors corresponds to different signals. 
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Figure 5; (Color online) Matrix plot represents values of window NRPDE for (a) NHP, (b) CHFP. x-axis represents window index 
and y-axis represents the samples. In each case, 1000 window is considered with length 10000. (c) mean NRPDE are drawn for 
NHP (blue) and CHEP (red). Each color in the matrix plot represents values of NRPD of each signals with the corresponding 
windows. 
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values of NRPDE corresponds to higher complexity in the dy¬ 
namics. So, window NRPDE infer that ECG signals of NHP 
groups always possess higher complexity in the dynamics. It 
is also observed that, window NRPDE for CHEP is less than 
the same for NHP (Pig.5b). This implies that, there exists less 
complex dynamics in the ECG signals of CHEP Thus it is im¬ 
portant to find a threshold to distinguish NHP and CHEP. The 
mean of window NRPDE is an effective tool in this context. 
In Pig.5c, we have shown the changes of mean window NR¬ 
PDE with respect to the window index. It is observed that, the 
values of mean window NRPDE in all cases of NHP lie above 
0.82 and thus whenever mean window NRPDE comes below 
0.82, the cardiac dynamics tending to be less complex in na¬ 
ture which may lead to congestive heart failure (Pig.5c). On 
the other hand, the values of mean window NRPDE is less than 
0.52 in all cases of CHEP. So whenever mean window NRPDE 
is found to be greater than 0.52 for some CHEP, the cardiac dy¬ 
namics becomes more complex in nature which corresponds to 
healthy cardiac dynamics. Prom the clinical perspective, this 
may indicate the recovery of the patient from congestive heart 
failure. It is also observed that the the difference between the 
mean window NRPDE of NHP and CHEP {d) maintains a fixed 
distance of 0.3 in all cases. All these results well correlate with 
clinical and experimental biomedical observations [2]. In this 
context, it may be noted that the similar observations for time 
series of consecutive heartbeat intervals for NHP and CHEP 
have also been observed in case of multiscale entropy (MSE) 
method [29,30]. However, the present method of window NR¬ 
PDE is obtained from the phase space of the given signal, while 
MSE is defined for the signal itself. As we know that the phase 
space gives the long term behaviour of the signal, the window 
NRPDE method is more robust than MSE that depends on the 
short term behaviour of the signal. Thus it is always expected 
that window NRPDE method will produce more unbiased result 
in the classification of NHP and CHEP 


4. Conclusions 

In this article we have studied the nature of complexity of 
ECG signals for normal healthy persons (NHP) and congestive 
heart failure patients. We have collected data from Physionet 
signal archive for 20 normal persons and 15 congestive heart 
failure patient (CHEP). The intention is to identify and distin¬ 
guish normal healthy persons and congestive heart failure pa¬ 
tients in terms of complexity. Complexity is recognized from 
RP of the corresponding phase spaces. In fact, probability den¬ 
sity of recurrent times, i.e; RPD have been calculated from RP, 
to observe an overall view of complexity of the dynamics. The 
order of complexity have been measured by NRPDE, the corre¬ 
sponding dynamics was verified well with known chaotic mod¬ 
els, stochastic systems and power noise. It has been observed 
that the dynamics of ECG for a NHP is more complex in na¬ 
ture compared to the same for a CHEP It is also observed that a 
congestive ECG follows a deterministic nature whereas the dy¬ 
namics of a healthy ECG is random. Einally a threshold value in 
the mean window NRPDE has been found out to distinguish the 


border line between NHP and CHEP. Though further investiga¬ 
tion into the physiological changes in ECG of NHP and CHEP 
is needed to consolidate the present study, our work reveals an 
alternative and at the same time a befitting procedure for the 
classifications of ECG signals of NHP and CHER Possible uses 
of this tool in clinical settings include the early detection of the 
life threatening congestive heart failure. 
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